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^^ ' Abstract. We introduce the Vainshtein mechanism which plays a crucial role in 

massive gravities, as well as in related theories such as Galileons and their extensions. 
This mechanism, also known as k-mouflage, allows to hide via non linear effects - 
typically for source distances smaller than a so-called Vainshtein radius which depends 
KjQ. on the source and on the theory considered - some degrees of freedom whose effects are 

then only left important at large distances, e.g. for cosmology. It is introduced here 
in non linear Fierz-Fauli theories (massive gravities), including the dRGT theories, in 
^ , their decoupling limits, as well as in other models such as DGF model or generalized 

Cp ' Galileons. This presentation is self-contained and before discussing the Vainshtein 

mechanism we introduce some useful results and concepts concerning massive gravity, 






j~^ , such as the vDVZ discontinuity, the decoupling limits or the Boulware-Deser ghost. 
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1. Introduction 

The idea to give a mass to the graviton is not new and has been investigated by many 
authors since the first years of General Relativity (for reviews see this volume and 
[U 12]). It recently regained popularity after the invention of the Dvali-Gabadadze- 
Porrati (DGP) model [3] and the discovery of its interesting cosmology [U [5]. Indeed, 
on the one hand, this model was the first where a large distance modification of gravity 
was shown to lead to cosmic acceleration even in the absence of a cosmological constant; 
on the other hand, the DGP model shares many properties which can be expected from 
a theory of massive gravity. The connection between the DGP model and massive 
gravity can of course be related to the fact that, in the DGP model, gravity is just 
mediated by a continuum of massive Kaluza-Klein gravitons due to the fiat and higher 
dimensional character of the bulk space-time. It is also amusing to note that there 
appears to be a historical connection between the very introduction of the cosmological 
constant (and hence the standard explanation of the later discovery of the accelerating 
expansion [6l [7]) and the idea to give a mass to the graviton (Einstein hoping that a 
cosmological constant would induced what is nowadays called a Yukawa decay of the 
gravitational potential j8]). 

The simplest theory for a non self-interacting massive graviton is known as the 
Fierz-Pauli theory [HI HO] . It suffers from a pathology known as the vDVZ discontinuity 
[TT] as will be introduced below. This is enough to rule out such a theory from basic solar 
system tests of gravity. However, soon after the discovery of the vDVZ discontinuity, 
a way out was suggested by Vainshtein [12], relying on a non linear extension of the 
Fierz-Pauli theory. This proposal of Vainshtein was however lacking a proper proof, as 
was underlined in particular in Ref. [13] which appeared soon after Vainshtein's paper in 
1972. The later reference also discovered another pathology of the non linear Fierz-Pauli 
theories, the so-called Boulware-Deser ghost, that we will also introduce below. The 
situation remained unchanged until the advent of the DGP model where the Vainshtein 
mechanism was re- introduced (the DGP model featuring also the vDVZ discontinuity), 
and where new arguments were given in favour of the validity of this mechanism [13] . 
Following this, a number of works tried to prove or disprove the Vainshtein mechanism in 
DGP or other simpler models related to massive gravity (see references in the following). 
The purpose of this work is to introduce with some detail this mechanism and discuss to 
which extend it can be considered established in the light of these recent developments, 
stressing that this mechanism plays a crucial role for phenomenological applications of 
massive gravity or of its close friends (such as Galileons, covariant Galileons and their 
generalizations [iSl [lEl [HI [E] ) . 

We first, in sections |2] and [3l introduce massive gravity theories: Fierz-Pauli theory, 
its simplest non linear generalizations and their generic properties. The results presented 
in these sections [2] and [3] are mostly known for more than 40 years (with the exception 
of the decoupling limit and associated strong coupling introduced in section 13. 3p and 
serve as a background for the following. In the way we present them, however, we 
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will use some more recent works which allowed to discuss them in a simple and clear 
way. We then, in section H] present more recent developments: the DGP model and 
the de Rham-Gabadadze-Tolley (dRGT) theory of massive gravity [191 EOl EI]. The 
latter theory was built specifically not to suffer from the presence of the Boulware- 
Deser ghost, a pathology that was once thought to be unavoidable in simple non linear 
extensions of Fierz-Pauli models (dRGT theory is presented in section |1]2]). In section |5l 
we introduce the Vainshtein mechanism and the recent progresses that have been made 
about its understanding in the various theories introduced before. 

2. Massive gravities, generalities 

2.1. Fierz-Pauli theory 

Fierz-Pauli theorjiji is the only consistent non self-interacting Lorentz invariant theory 
for a massive spin 2 P [10]. It can be defined by the following action defined on fiat 
space-time with canonical metric g^i, = rj^^, 



SpF = Ml I d^x 



\ {d,K,f + \ {d,hf - \ id,h) {d^W,) 



+\ {d^Kp) {d'h^n - \m' {h,.h^^ - h') + MfT^^h^^'' 



(1) 



where h^y is some rank-2 covariant tensor, m and Mp are mass parameters, the indices of 
h^y are moved up and down with the metric rj^y and h = h^yrj^'-', while hfj^^ is coupled to 
some energy momentum tensor T^^,. The graviton hf^y has a mass term in the Fierz-Pauli 
form, that one reads from the above equation as given by 

SpF,m = - ^M^m^ f d^x [h^yh^^- - h^) . (2) 

Any other combination of h^yh^'^ and /i^ would lead to instabilities [9l[T0]. Note also that 
this mass term explicitly breaks gauge invariance. The first four terms of action ([1]) are 
simply obtained by expanding at quadratic order into h^y the Einstein-Hilbert action 
around fiat space-time. Varying action ([1]) with respect to h^y, we get the following 
equation of motion 

£p,y = -9"^^ ^^^"' " ^^^"""^ ^ Mp'^Tf.y, (3) 

where £fj,y is the linearization around ri^^y of the Einstein tensor G^y and is given by 

= - \d^dyh - ]ph^y + ]^dpdX + \dpdvh'^. - IvpAd'd'^h,^ - Dh) 

I In this work, we will keep the name " Fierz-Pauli theory" to describe the quadratic theory introduced 
in this subsection. This theory has linear field equations. By contrast, we will call "Non-Linear Fierz- 
Pauli theory" - henceforth NLFP - a non linear completion of the later theory of the type introduced 
in the following section [31 
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It is easily checked that, as a consequence of Bianchi identities, £^^ is divergenceless. 
Hence, if one further assumes that T^^, is conserved with respect to the background 
derivative d, i.e. that one has d^T^^, = 0, we get (taking the divergence of the field 
equations ([3])) for non vanishing graviton mass m 

d'h,^ = d^h. (4) 

This represents 4 first order equations which eliminates as much degrees of freedom out 
of the 10 a priori independent components of h^^. Taking one more derivative of the 
above equation, we have 

d'^d^'h^^ -nh = 0. (5) 

The left hand side of the last equation is nothing but the linearized Ricci scalar. Eq. 
([5]) can be fruitfully used in the trace (with respect to rj^^) of the field equations ([3]), to 
get 

Hence the trace h of h^v is determined algebraically and does not propagate. We will 
use this equation below, but just notice at this point that, if one considers the theory in 
vacuum, equations (jlj) and ([6]) imply that a massive Fierz-Pauli graviton in vacuum is 
transverse and traceless, and hence contains 5 degrees of freedom (the same conclusion 
can be reached by a rigorous hamiltonian counting, as we will recall later). 

2.2. The vDVZ discontinuity 

Plugging (jlj) and (Q into Eq. ([3]) we get easily that 

Using then a Fourier decomposition of /i^,^ (and using similar notations for T^j,) 

V(xO = T^ J d^k e^'^"^" V ik") , (8) 

we get from ([7]) the expression of the propagator (in Fourier space) D)^pa , such that 

rppa 

h =D^^^^)1 (Q) 



pupa ^2' 



given by 



^pupa = ^2 _|_ ^2 \VppVcru + VpuVcrp — -^VpaVpu ~ ^'7po- ^3 ) ' V^^ 

where k"^ = k^k^^. This should be compared to the propagator of a massless graviton 
Dpjpa which reads in Fourier space 

^t<^7^^ ~ p (^PM^^^ + VpuVap - rjpaVpu) , (11) 

where the symbol ~ means that gauge dependent and momentum dependent terms 
are omitted. Notice the crucial difference between expressions (llUp and (II ip lying in 
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the third term in the parentheses on the right hand side of these expressions. This 
difference in the coefficient in front of rip^rj^i, is independent of the mass of the graviton 
and is at the root of the so-called vDVZ discontinuity (vDVZ standing for "van Dam- 
Veltman-Zakharov" ) which states that, roughly speaking, however small the graviton 
mass, Fierz-Pauli theory leads to different physical predictions (such as light bending) 
from those of linearized General Relativity [11]. There are various ways to see this. Let 
us here (following e.g. [22j) consider the formal tree level amplitude A between two 
conserved current T^u and Sf^,^, defined as 

A=M], fd^x 5^"(x)V (T) (a;), (12) 

where hf^u (T) is the tree level graviton field generated by the conserved source Tp^j, and 
given by 

h^,{T){x) = Mf j d'x' D^,p, {x - x') T"'^ (x') , (13) 

{Dfj,upa being the massless or massive propagator). The amplitude A is easily obtained 
in Fourier space, using the expressions ( TTOj) and ( ITTj) for the propagators. One gets 
respectively in the massless and massive cases, and in the large k limit {k ^ m) 

A^'^='^ = fd'k^ (s>'%. - ^Sf\ (14) 

A^"^^'^ = fd'k^ (s'^'Tp, - ^Sf] . (15) 

Considering then non relativistic currents - such that T^ oc diag(Mi, 0, 0, 0) and 
Sj^ oc diag(M2, 0, 0, 0) - separated by a distance small with respect to the graviton 
Compton wavelength (which diverges as m goes to zero), the amplitude due to the 
exchange of a massive graviton is given approximately by 

^i^m^U-«>^UAMiEi, (16, 

so that the massive amplitude stays different from the massless one, however small 
the graviton mass. This translates into a similar discrepancy in the non relativistic 
potentials between the two same sources, the potential of the massive theory being 
larger by a factor 4/3 and it reflects an extra attraction in the massive theory with 
respect to the massless theory. This extra attraction can be attributed to the exchange 
of the helicity zero polarization of the massive graviton (which, as we recalled above, 
has 3 more polarizations than the massless one). It can be eliminated by redefining the 
Newton constant of the massive theory with respect to the massless one, assuming, e.g., 
that one measures the Newton constant by some Cavendish experiment (indeed, if one 
does not do such a rescaling, the Newton constant of the massive theory would be given 
by 4/3 X Mp^/lGvr = Mp^/127r). However, with such a rescaling, the discontinuity will 
then reappear in other observables, like the light bending. The latter will then be 25% 
smaller in the massive case than in the massless one [llj, which is much too large to be 
compatible with current measurements of the light bending by the sun [231 El]- 
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As such, the vDVZ discontinuity is enough to rule out from standard solar system 
tests of gravity any theory where it appear4§|. It must be cured if one wants to make 
such a theory viable. The Vainshtein mechanism (first introduced in Ref. [12]) allows 
in principle to get rid of the vDVZ discontinuity, as will be explained in the following. 
It relies on non linearities in the field equations which are absent (by definition) in 
the Fierz-Pauli theory. It is also clear that a gravity theory able to approach General 
Relativity in high curvature regimes must be non linear, independently of the Vainshtein 
mechanism, and they are then several good reasons to consider non linear completion 
of the Fierz-Pauli theory. 

Historically, the non linear completion that played the most important role for the 
recent developments about massive gravity, is certainly the Dvali-Gabadadze-Porrati 
model (henceforth DGP model) [3]. This model will be briefly discussed below in section 
14.11 and, as we said, it is this model and its cosmological consequences, that, in the years 
2000, lead to a very strong renewal of interest about massive gravity, related theories 
(such as, e.g., Galileons [15]) and the Vainshtein mechanism (as discussed in [H]). 
However, the DGP model is far from being the simplest non linear extension of Fierz- 
Pauli model one can consider, since in particular it contains a non countable infinity of 
massive gravitons. We introduce in the next section a family of much simpler non linear 
completions, that we will call here and henceforth "Non Linear Fierz-Pauli" theories 
or NLFP. The Vainshtein mechanism was originaly introduced in a theory close to one 
of this familjjljj and the recently discussed massive gravity of de Rham, Gabadadze and 
Tolley [T9l l20l [21] also belongs to this family. Note that the NLFP theories were first 
discussed in relation with strong interactions [31] . 

3. Non Linear Fierz Pauli gravity 

3.1. Action and equations of motion 

An obvious way to generalize in a non linear way the Fierz-Pauli theory considered in 
the previous section, as well as to stay close to General Relativity (henceforth GR), is 
to consider a theory with a dynamical metric f^^j, and the same kinetic term as the one 
of GR, hence given by the usual Einstein-Hilbert action reading 

Seh = ^ (/ d'^V^ (R - 2A)) , (17) 

with a possibly non vanishing cosmological constant A, and to add to this term a mass 
term. This mass term should be such that when expanded at quadratic order around 
some suitable background for g^^, namely that given by a flat metric rj^^^, it reduces to 

§ Note however that the discontinuity does not appear if the background is maximally symmetric with 
a non vanishing curvature [331 HSl US] and that this result can be extended to more general situations 
[211 EHl ES] ■ However, depending on the value of the mass of the graviton, the theory can be non unitary 
when background curvature does not vanish [25 . 

II It was pointed out in [30 that the theory used by Vainshtein does not fall strictly speaking into this 
class, but it does not matter for the discussions of this paper. 
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the Fierz-Pauli form ([2]) where h^^, is now viewed as the dynamical field representing 
the fiuctuation of the metric g^^, around its background expression. Ideally, the mass 
term for g^i, should only depend on g^i, itself in a non derivative way. However, the only 
such non trivial term corresponds to a Lagrangian density proportional to the volume 
element \/—g, hence to a cosmo logical constant. Moreover, such a term, if expanded 
around some arbitrary background metric contains a tadpole (in particular when the 
background metric is the Minkowski metric), and does not have the Fierz-Pauli form 
at quadratic order. In fact, it is clear that it does not give any mass to the graviton 
since this term does not break general covariance and hence leads to a theory with two 
propagating degrees of freedom. 

Hence, the sought for mass term requires the introduction of some extra field besides 
the metric (7^1,. One possibility we shall consider here is to introduce an extra metric 
ffj^iy that will be taken to be fiat and non dynamical. As will be explained below, such 
a theory can equivalently be seen as that of a dynamical metric g^y and four scalar 
fields. The two metrics f^^p and g^v will be taken to be non derivatively coupled via an 
interaction term Sint[fig]- This term will be chosen such that (i) the theory is general 
covariant under diffeomorphisms (common to the two metrics), (ii) it has fiat space-time 
as a solution of the field equations for g^^, and (iii) that when one expands g^y to second 
order around the canonical Minkowski metric rj^y as g^i, = tj^i, + h^^ and let f^^ to have 
the canonical Minkowski form r]^^, the potential at quadratic order for h^y takes the 
Fierz-Pauli form ([2]) . There is much freedom in the choice of such an interaction term. 
For example, the following two possibilities have be considered respectively by Boulware 
and Deser in Ref. [13] and by Arkani-Hamed et al. in Ref. [32] 

5^ ^ - ^m'Ml I d'x v^ H,,H^, {r T - FT) (18) 

S^ ^ - ^m'MJ, J d'x v/^ H,,H^, {g^^g''^ - g^'^g^^) , (19) 

where f^'' and g^" denotes respectively the inverse of the metric f^y and g^u, and H^u 
is defined by iJ^jy = g^^ — f^u- More generally, one can consider theories where the 
interaction term Sint is not of the above forms (llStllQp l^. but is chosen such that it obeys 
properties (i), (ii) and (iii) above. We will later introduce the interaction terms specific 
to the dRGT theories [T^l [201 [2T] which also differ from the ones above. Note that since 
we have two metrics at hand, there is some ambiguity on how to move indices up and 
down. Here, when necessary, indices will be moved with the dynamical metric g except 
for indices of the other metric f^y itself. I.e. f^^ is defined as the inverse of the metric 
f^y and hence it is it not given by g^^g'^^fap (this means in particular that H'^'^ is not 
given by g^" — f^^). Following the notations of Damour et al. [33], one notes that the 
interaction terms considered so far all have the form 

41 = -\m'Mlj d'xV''^\g,f) ^ -\m'Ml j d'x./^gV^-\g-H) (20) 
^ We keep the numbering of the interaction terms (|18m9p of Ref. [53] . 
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with V^"-* ((?,/) = y^—g V^°'\g~'^i) a suitable "potential" density associated with the 
scalar function V^^\ In fact one can show that this form follows necessarily from the 
assumption that the metric / and g interact in a non derivative way in the interaction 
term Sint [/, g] , and that the theory is invariant under diffeomorphism ^33j . This later 
invariance is represented as usual as the transformations acting on the metrics of the 
following form 

g^^u{x) = d^x"'{x)dy^{x)g'^^{x'{x)) , 

UA^) = d,x"^{x)d,x'^ix)f:,A^'i^)) , 

and under which the quantity V^"''^ transforms as a scalar. Introducing matter, there 
is also much freedom as far as choosing the metric to which matter couples. Indeed, 
e.g., an infinite family of inequivalent metric can be built from the two metric at hand 
f^u and g^u, and one could decide to couple matter minimally to one arbitrary metric 
in this family. Here, we will consider the simplest case where matter is assumed to be 
minimally coupled to the metric g, hence the total action of the theory we will consider 
here is given by 

S = j d'x^^g (^R + L)\+Srnt[f,gl (22) 

where in the above action, Lg denotes a generic matter Lagrangian with a minimal 
coupling to g^j_v (and not to the metric f^v)-, and we have included a possibly non 
vanishing cosmological constant in the interaction term Sint[f,g\- 
The field equations, derived from action ( l22l) . read 

MJ,G,,= {T,, + T^,), (23) 

where G^^ denotes the Einstein tensor computed with the metric g, T^^ is the energy 
momentum tensor of matter fields, and T^^ is the effective energy momentum tensor 
coming from the variation with respect to the metric g of the interaction term Sint- It 
depends non derivatively on both metrics / and g and is defined as usual as 

TU^) = - At. L A nt[f,9]- (24) 

\'~gog^[x) 

It is then easy to check with these expressions that the equations of motion (123!) reduce 
indeed to the Fierz-Pauli equations ([3]) at linearized level. A simple, but non trivial, 
consequence of equations fl23|) is obtained by taking a g-covariant derivative V of both 
sides of the equations; one gets, using the Bianchi identities and the conservation of the 
matter energy momentum tensor, the constraint 

V^T^t = (25) 

which the effective energy momentum tensor should obey. 

We will call a theory having an action of the form fl22|) discussed above (where 
we recall that the interaction term obeys conditions (i), (ii) and (iii) above) a Non 
Linear Fierz-Pauli (NLFF) theory. We will now introduce a pathology, first discussed 
by Boulware and Deser jT3], and once thought to be present in any NLFP theory. 
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3.2. The Boulw are- Deser ghost 

It is here convenient to contrast a Hamiltonian analysis of the Fierz-Pauh theory with 
one of a generic NLFP theory, the starting point being a 3 + 1 decomposition of the 
graviton or of the metric. 

Since the kinetic term of the hnear Fierz-Pauh theory is the same as that of 
hnearized General Relativity, one can use first well known results concerning GR. 
In particular, one sees that neither /iqo nor h^i (with i = 1,2,3 spatial indices) are 
dynamical degrees of freedom since their canonical momentum vanish identically in 
massive gravity as well as for a massless graviton and /iqo and hoi are Lagrange multiphers 
in the kinetic term obtained by expanding the Einstein-Hilbert action. In the mass term, 
however, /ioo and /iqj play quite different roles. Indeed, the mass term 02]) reads 

- \mW I d'- {h,K, - 2h^h,. - K^„ + 2/,,/,„,} , (26) 

and it hence appears that that /iqo is a Lagrange multiplier for the entire action since 
it also enters linearly in the mass term. As a consequence, the field equation for /iqo 
generates a constraint which reads 

V^/iii - /iij,y oc m^hii, (27) 

and allows to eliminate one extra degree of freedonO. In contrast, the field equations 
for the /loi, which appear quadratically in the mass term, do not eliminate any degree of 
freedom (in contrast to what happens in the massless case), they determine the /iqj in 
terms of the other dynamical components. Together, this leaves a total of 5 propagating 
degrees of freedom, in agreement with the discussion given above in section 12. 1[ 

Let us now see how those results are modified in the case of a NLFP theory. In the 
massless case, that is to say for General Relativity formulated d la ADM [31] , Lagrange 
multipliers associated with diffeomorphism invariance are the "lapse" iV and "shifts" 



A^*, respectively defined as A^ = 1/y — (7"° and Aj = (^oi in terms of the components 
of the metric g^^. They generate (first class) constraints which eliminate 4 out of 
the 6 possible dynamical degrees of freedom, leaving the well-known 2 polarizations of 
a massless graviton. The addition of a NLFP "mass term" such as f fT8]) - flT9]) modifies 
however notably the nature of A^ and A^* (as it breaks invariance under diffeomorphisms) . 
For example, the action with the mass term (fT8l) reads in the first order formalism (after 
a convenient renormalisation of rra^) 

Ml I d^x { [n'^gij - NR° - NiR') 

-m" {hijhij - 2NiNi - hiihjj + 2hii {1~N'' + Nkg'^'Ni) ) } , (28) 

where the vr*-' are conjugate momenta to the gij and R^ and i?* are respectively the 
Hamiltonian and momentum constraints of GR (generated by the lapse and shifts) and 

+ To state things more rigorously, on can show that this constraint is second class and generates an 
extra constraint of the same class. These constraints together kill two Hamiltonian degrees of freedom, 
i.e. one Lagrangian degree of freedom. 
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cubic terms simply follow from the standard definitions of A^ and A^j recalled above. As 
Boulware and Deser first pointed out [13j neither A^j nor A^ are Lagrange multipliers 
of the non linear theory. Hence, the number of propagating degrees of freedom is 
generically 6 and not 5. Boulware and Deser also argued that the reduced Hamiltonian 
for the 6 physical degrees of freedom is in general unbounded from below, and this can 
indeed easily be checked explicitly in some cases such as that of action f l28|) (see e.g. 
|35j). Given the unboundedness-from-below nature of the Hamiltonian the extra-mode 
is usually called the " Boulware-Deser ghost". Note also that (as discussed e.g. in [35] ) 
one can also understand the presence of the Boulware-Deser ghost in a covariant way 
using the observation that the constraint (E]) is lost in generic NLFP theories and is 
replaced (at cubic order), schematically, by a quadratic equation of the form 

Dh^ + m^h oc Mp^T (29) 

(while a constraint such as (j4j) still holds in the form of (l25l)). Note that this discussion 
is modified for dRGT theories which were precisely built to eliminate the BD ghost. 

3.3. Strong coupling and decoupling limit 

It was first noticed in [14j, in the context of the DGP model, that the Vainshtein 
mechanism can only work at the price of a strong coupling in the considered theories. 
A subsequent work [32j introduced a powerful method to extract from a generic NLFP 
theory a simple theory which captures many of the crucial properties concerning this 
strong coupling and its relation with the Vainshtein mechanism. The purpose of this 
subsection is to introduce this method and some of its outcomes. 

A good starting point is to notice that the gauge invariance f l2T|) can be used to 
write the background flat metric / in various coordinate systems. Starting from a 
given gauge, with coordinate X^, and the / metric in the form of /ab(X), it might 
be desirable to change the gauge, but keep the change of coordinate explicit in the / 
metric. Namely, the action considered takes the form of f l22|) . but with f^u{x) now given 
by the expression 

/^.(x) = d^X\x)da''{x)fAB {X{x)) , (30) 

while g is kept as g^u{x). The quantities X^, which then appear explictly in the action 
of the theory, can be considered as a set of four new dynamical scalar fields, which 
are analogous to the Stuckelberg field used to restore gauge invariance in the Proca 
Lagrangian |32| |36]. 

With this in mind, the initial gauge, where g and / assume the form gAB and fAs 
is usually called a "unitary gauge", i.e. one where the Stuckelberg fields X"^ are gauged 
away. Note that the metric f^u in a non unitary gauge can also be thought as the 
puUback, via the "link field" X"^(x), on the space-time manifold 777.4, with coordinates 
x^., of the metric Jab living in an other abstract manifold A^4 with coordinates X^ [32] . 
Usually, the unitary gauge is chosen such that in this gauge the metric /as takes the 
canonical Minkowski form tjab = diag(— 1, 1, 1, 1). In the non unitary gauge, the action 
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is one for a theory with g^j^^y and X^ as dynamical fields. Obviously, the equations 
of motion for g'^'^ lead to the same equations as in fl23l) where /^j^ is given the form fl30|) . 
On the other hand, it is not difficult to show that the field equations for the X^ are 
equivalent to the Bianchi identities ( 125|) provided that the mapping X'^(x) is invertible 
(for an exphcit proof see e.g. [37]). 

The authors of Ref. [32] have further developed the above mentioned analogy 
between X^ and Stuckelberg fields doing a " Goldstone boson" expansion of the action 
( l22l) around a unitary gauge. Considering some background solution for (yf^y (defined as 
g^j^) and X^{x) defined (the metric /ab being kept fixed) as 

X^ix) ^ Sy, (31) 

Ref. [32] introduces the "pion" fields n^ as 

X^(x)=Xo^(x) + 7r^(x), (32) 

and further does a " scalar- vector" decomposition of the tc^ in the forno 

7r^(a;) = 6^ (A^ix) + r^^^S.^) . (33) 

If one inserts this decomposition into action ( l22l) . and expands around fiat space-time 
writing g^^^, = rj^y + h^y, we obtain an action for the dynamical fields h^^{x), A^{x) and 
4>{x). Since A'^{x) and 0(x) only enter in the metric f^^ (via expression ( l30l) ) the only 
term in action (|22|) which depends on A^{x) and 0(a;) is the interaction term Sint[f,g], 
and one has (where no term has been neglected in the expression below) 

H^i, = h^y - d^Ay - dyA^ - 2d^dy(j) - d^A^dyA" 
-d^d^<pdyd''<P-dyA^d^d^<f)-d^A''dyd,<f). 

Following again p2] , we can obtain canonically normalized fields 0, A and /i^,y by defining 

h^y = Mph^y, A^ = MpmA^", (p = Mprn^cj). (35) 

Inserting (IMl) with (l35ll into Sint[f,g\, keeping the lowest order in /i^,^(x), A^^{x) and 
0(x), and using the redefinition 

hf^y = h^y - r]f,y4), (36) 

one obtains the following quadratic action 

S=- f d'^xhh^"'d^dyh - 2h'"'dyd^h''^ + h^^Uh^.y - hUh 

+ m^ {h'' - h^yh^'"^ - F^yF'^'' - Am (hdA - h^yd^A"^ (37) 

+ 60(n + 2m2)0 - m^hcj) + 2m0(9i| + -^h^"' - i— 

J 2 ikf p 2 ikf p 



where h = h^yfj^^^ , dA = c^^jA^, T^y is the matter stress-energy tensor, T = T'^^rj^y, and 
indices are moved up and down with the metric rj^y. The peculiarity of the above action 

* This form is in fact not explicitely the one given in [32| but it seems to be what is done there 
imphcitely. See j37j for a discussion about the associated subtleties which do in fact matter when one 
discusses the terms beyond quadratic order. 



An introduction to the Vainshtein mechanism 12 

is that while A^^ acquired directly a standard kinetic term, did only via a mixing 
with h^^ (that was demixed through definition fl36|) ) [32], this being entirely due to the 
structure of the Pauli-Fierz mass term (jj]). The left over cross terms between 0, A and 
h^i, can be cancelled by adding an appropriate gauge fixing to the action (see [38]). 

Expanding the action in 0, A and /i^i, to next orders, one sees that has in general 
cubic self interactions suppressed by the energy scale 

A5= {m^MpY^' . (38) 

When these interactions are presenlij, they are the strongest interactions among the 
fields 0, A and hf^^, in the limit where m ^ Mp. One can take a decoupling limit 
(henceforth DL) defined as 

Mp — i- oo, m — )• 0, As ~ constant, T^^/Mp ~ constant, (39) 

in order to isolate this interaction. In this limit, the action one is left with for is of 
the form 

S = \j d'x { ^000 + -1 [a (00)3 + /3 (n0 0,^, ^'n] - -^/~A ' (40) 

where a and /3 are numerical coefficients that depend on the interaction term Sm/lf, g]iS 
For example, the potential flTSj) leads to a = — /3 = —1/2, while the potential flT9|) leads 
to the opposite case a = —(3 = 1/2. In contrast, in the DL, the other fields h^,^ and A^ 
are just free. 

The equation of motion deriving from action (HOj) is 



A5 



3« D (00)' + /3 D (0,^, c^''^'^) + 2/3 5^9, (n0 0'^^ 



^ r.(41) 



As will be discussed in section [5.1. li the Vainshtein mechanism can be easily read off this 
equation. Moreover, the fact that this equation is of fourth order signals that action PUI) 
propagates in fact two scalar modes, one being ghost like. This ghost can be interpreted 
as the Boulware-Deser ghost and the DL provides a powerful tool to investigate the 
presence of this ghost in a given theory (as first argued in Refs. |2Sl 1^)- 

4. From DGP to dRGT gravity 

4.I. DGP gravity in brief 

The simplest DGP model [3j is a five dimensional (5D in the following) brane-world 
model, and as such, it describes our four dimensional (4D in the following) Universe 
as a surface embedded into a 5D bulk space-time. Standard matter fields are thought 
as being localized on this surface, while the gravitational fields are living in the whole 

jj An appropriate choice of the interaction term Smf [/, g] can remove cubic (and some others) self 
interactions of cj) [32] . 

ft Note that in general, the cubic term for (p is given by some linear combination of the three terms 
(D^)'^, 0(f) (f)_^i, (f)''^'^ and ^^^i, 0''"^(^'J^, but an integration by part can always be used to reduce the 
number of independent term to two, as shown in Eq. (j40| . 
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bulk space-time. The characteristic feature of the DGP model lies in the gravitational 
dynamics. The latter is first governed by a bulk gravitational action, which is the usual 
action for 5D gravity, namely 

S{5) = ^J d'X,/gl^R^,^, (42) 

(5) 

where i?(5) is the 5D Ricci scalar computed from the 5D metric (7^^ and M(5) is the 
reduced 5D Planck mass. To account for the brane, one adds to this action a term of 
the form 

3(4) = J d^x,/g^ (^R^^^ + £(M)) , (43) 

where £(m) is a Lagrangian for brane localized matter (that is to say baryonic matter, 
dark matter, ...), and i?^'^) is the Ricci scalar of the induced metric gliv on the brane. 
It is this term, depending on R'^^\ that is responsible for all the peculiarities of the 
gravitational phenomenology of DGP gravityljj The induced metric g\jiv is the metric 
experienced by the matter we are made of, it is defined by giiy = d^X'^dyX^g^^ , where 
X°'{x^) are defining the brane position in the bulk (X" being bulk coordinates and x^ 
coordinates along the brane world- vo lume) . 

In this model, the gravitational potential between two static sources, separated by 
a distance r, interpolates between a 4D 1/r behavior at small distances and a 5D 1/r^ 
behavior at large distances, as shown in reference p5j. The crossover distance r^ between 
the two regimes being is given by 

However, gravity does not reduce to a Newtonian potential, and in the DGP model, 
from a 4D point of view, gravity is mediated by a continuum of massive gravitons 
(so-called Kaluza-Klein modes), with no normalizable massless graviton entering into 
the spectrum. As a consequence, the tensorial structure of the graviton propagator was 
shown to be the one of a massive graviton and the model shares some properties with the 
previously introduced NLFP theories. As we said, DGP gravity features in particular 
the vDVZ discontinuity. 

A decoupling limit can be obtained in DGP gravity as well. In the scalar (i.e. 
helicity 0) sector, the obtained theory has a Lagrangian reading [10| W\\ 

C,,, = 3000 - -^{d,~4>fU4^ + ^, (45) 

where the strong coupling scale is given by Kdgp = (Mp/rlY^^. Again the presence of 
the cubic operator in the above Lagrangian allows to understand easily the Vainshtein 
mechanism and associated scalings. Note however that the field equations deriving from 
this Lagrangian are just second order in line with the expectation that the DGP model 

f This term is expected to be present in generic brane world constructions. It is the hierarchy of scales 
between Mp and Mf^-^ that is the real distinctive feature of DGP gravity. 
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does not contain any Boulware-Deser ghost [55]. This remark is also at the root of the 
introduction of the Gahleon family [T5] . 



4-2. dGRT gravity in brief 

After the discovery of the BD ghost, the question arose if this ghost could be eliminated 
by a suitable choice of NLFP theory. This question was raised in the article of Boulware 
and Deser of 1972 [T3], and later, using the decoupling limit, in Ref. [39] which concluded 
that it was in fact not possible. The same question has been recently reexamined by de 
Rham, Gabadadze et Tolley [191 1201 121], who, using first the decoupling limit, built a 
family of theory which appeared (in the DL) to be devoid of the Boulware-Deser ghost. 
This was later confirmed by a Hamiltonian analysis of the untruncated theory [^ HS] 
and several other works [S], even though some authors still disagree [IHl Il6l Wf\ . 

dRGT theories are NLFP theories for which the interaction term S'j„t[/, g\ of action 
( l22l) has a special form. Using a parametrization given in references [42|l48]. we introduce 
functions e^ defined for an arbitrary n^iu matrix[j|X^j, represent elementary symmetric 
functions of the eigenvalues of X, and given, for A; = to 4, by (cf. e.g. [12]) 

eo (X) = 1, ei (X) = [X], e, (X) = \ {[Xf - [X^]) 

63 (X) = 1 {[Xf - 3[X][X2] + 2[X3]) (46) 

64 (X) = 1 ([X]^ - 6[X]2[X^] + 3[XT + 8[X][X3] - 6[X1) 

where [X] designate the trace X\ of the matrix X. For an arbitrary integer fc, one 
defines Ck by 

e,(X) = ix^^[,,...X^'=4], (47) 

where brackets [ ] around indices indicate the unnormalized antisymmetric sum over 
permutations. For a matrix X, Cayley-Hamilton theorem implies that 

det(X) = e„(X). (48) 

dRGT theories [HJ [201 [21] can be defined by an action of the form [42] 

k=A 
fc=0 

where Pn are arbitrary parameters, and the square root above is a matrix square root 
of the tensor g~^fl§|. One can notice, on the one hand that /3o describes a cosmological 
constant (which does not give any mass to the graviton), on the other hand that the 

I With /, a line index belonging to {!,..., n}, and J, a column index belonging to {l,...,n}, and n 
having so far no relation with the space-time dimension D. 

§ Note that this square root does not always exist. However, a vierbein formulation (such as given 
in [49]) of the theory allows to show that the existence of the square root does follow from the field 
equations, at least in a subset of theories. It also provides a simple way to count the number of degrees 
of freedom El [50]. 



S = Mljd^x./^ i?-m2^/3fcefc(v/^^ 
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term proportional to (^4 does not give any contribution to the field equations of g^^, 
since ^/^ 64 f Vy^j = \^~9 det ( \/g^^f) = \/~f- Hence, for D = 4 dimensions, 
we get a three parameter family of massive theories parametrized by Pk, with k = 1,2,3 
(which becomes a two parameter family once the mass of the graviton is fixed). It is 
easy to extend this construction to D dimensions by considering actions of the form 

k=D-l 



S = M^-2 I d^xy/^ 

k=l 



(50) 



where e^ are defined as in fj47|) . The interaction term Sint[f,g] of theories fH9|) - fl50|) is 
indeed of the form f l20|) . 

4.3. Decoupling limit of dRGT model 

A decoupling limit can also be obtained in the dRGT theory, as we now recall. It will be 
more convenient for that purpose to use an alternative form for the action ( 149|) . written 
in terms of the matrix /C, 



/C = I-Vg-if. (51) 

The dRGT action then reads, 

S = Ml I d'^x^^ [R + 2m^ (62 (/C) + 0363 (/C) + 0464 (/C))] , (52) 

with the following identifications, 

/3o = —12 — 80:3 — 20:4, /3i = 6 + 603 + 0^4, 
132 = -2- 4a;3 - 20:4, fi^ = 2(03 + 04) 

As in the case of NLFP, we introduce the "Goldstone boson" expansion, defined by ( 15U]) . 
f p2|) and f p3|l . The matrix /C defining the interaction term has then the form, 

/c;; = (5^ - ^/W^^thZ, (53) 

where H^^ is given by fl34|) . Expansion of f l52|) in powers oi H = H^ = {g~^H) gives, 
62 (/C) + ^363 (/C) + ^464 (/C) = ^ ([H]' - [H']) 



+ Y^{y[i^r + 2(^3 - «4)[i^F[i^'] + (1 - 2^3 + a4)[i^T 
+ 4 (1 - «3 + '^) [H][H'] - (5 - 4^3 + 2^4) [H']} + OiH'). 



Up to redefinition of constants, fl54j) coincides with equations obtained in [20] and 
[5T] . Substituting fl5^ with i^ expressed as (IMl) into fl52|) . one obtains the action 
as a series expansion in h^iy, A^ and 0. Because of the particular form of 62 (A^) in 
( 152|) . chosen to recover the Fierz-Pauli theory, the quadratic part of the action leads to 
the same expression, as in general NLFP, Eq. ( 1371) (up to an overall factor 2, due to 
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different normalizations in the actions). As we already mentioned before, it is possible 
to arrange the interaction term of a generic NLFP to remove the leading interactions 
~ [dipY /{Mpm^). In fact, the dRGT theory does exactly this: the interactions 
~ [d^Y / [M piTL^) cancel and the leading interaction is then suppressed by a higher 
scalqjjl 

A3= (m^Mp)'/'. (55) 

Therefore for the dRGT model we define the decoupling limit as 

Mp — 7- oo, m — 7- 0, A3 ~ const, T^^/Mp ~ const. (56) 

Since the vector mode does not couple to a source, we can set it consistently to zero, 
A^ = 0. Then in the decoupling limit (l56ll we obtain. 






JL 

where the following notations are introduced: (5 = 1 + 0:3, /3 = 0:3 + 0:4, 



(57) 






{51 



XS) = - ^6/^^e/%<l>„^$,., (59) 

^S^ = le,''''ef-'^ap<^p.<^,5, (60) 

with $^;^ = (j)^fj,i,, and e the totally antisymmetric Levi-Civita tensor. Note that as it 
was argued in [201 HH] the higher order terms in the decoupling limit are equal to zero, 
therefore we stopped the expansion in (154]) at the fourth order. The equations of motion 
obtained from 05 7p read, 

'^x(2) - Ax(3) = la 

q M^ A6 '"' Mp 

(62) 

Notice that the helicity-0 and helicity-2 modes are mixed at all orders. It is possible, 
however, to diagonalize the linear and quadratic mixings. Indeed, making here the 
nonlinear field redefinition in 057p |2U] (as opposed to the linear redefinition 0361) used 
before). 



£"^/i R - x(^) - — x(2) - --x(3) = -^ r6i) 



h^u = h^y - Tj^^cp - a — -3 — , (63) 



II Note that the cubic vector-scalar terms ~ dAd^(f>d^(f> are suppressed by A4 ^ A3, however, these 
terms cancel out up to a total derivative |52) . 
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and using the relations S^j^ {rjai3(j)) = —X^J and £^^ {da4>di3(j)) = —X^J one obtains up 
to total derivatives, 

S = I d'x[ - h^'S'^.^Kp + ^000 - ^0'>'"X« 



+ s;(^-''"-^^-^^-^''^^^)}- 



One can recognize the Galileon [151 US] kinetic terms in flMj) . Notice that the kinetic 
mixing does not disappear in the general case /3 7^ 0. Only in the special case /? = 0, the 
helicity-0 and helicity-2 modes are completely decoupled. It is also worth to note that 
couples directly to matter source in a "disformal" way. Varying ( IMj) with respect to 
h^y and respectively, we obtain, 

^.t V - ^4t^ = % (65) 

300 + l^'^^^l^ + ^ (f + f ) ^"^- + ^^^^4t^ 

5. The Vainshtein mechanism 

5.1. The Vainshtein mechanism in the decoupling limits and k-mouflage 

The Vainshtein mechanism was introduced by Vainshtein by looking at static spherically 
symmetric ansatze in a theory close to one belonging to the NLFP family discussed 
before. Such a theory, at large distances from a source, can be described linearly, so 
that all modes, including the ones additional to the two polarizations of a massive 
graviton propagate, resulting in large deviations from GR responsible for the vDVZ 
discontinuity. Vainshtein pointed out that, at some distance to the source, called today 
the Vainshtein radius, ry, the linear regime breaks down, and for r below ry the theory 
enters a nonlinear regime. The Vainshtein radius is a composite scale made out of the 
Planck mass, a scale specific to the theory considered (e.g. the graviton mass) and the 
mass of the source. Vainshtein also proposed to find the solution below the Vainshtein 
radius using an expansion in terms of the mass of the graviton, instead of that in terms of 
the Newton's constant that is used at large distance. However he did not show that the 
two expansions could be matched in an existing solution. That such a solution exists 
was only proven recently in Refs. [371 [531 [5l] . However, the Vainshtein mechanism 
was re- introduced before, in the framework of DGP gravity |T3]. There, indeed, exact 
solutions - however of cosmological type - were found jU [55] featuring for the first time 
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an interpolation between two regimes: one exhibiting properties similar to those of the 
linearized regime of NLFP theories and another recovering GRiill. 

The Vainshtein mechanism was then used and believed (and in some cases has been 
shown) to allow the recovery of General Relativity around massive bodies, and more 
generally in cases where the space-time curvature is larger than some critical value, for 
a wide class of non-GR theories of gravity. Roughly speaking, this is achieved by hiding 
extra degree (s) of freedom by strong kinetic self-coupling, so that they almost do not 
propagate. That it can work can already be seen in the decoupling limit, described in the 
previous sections, which allows to separate the different helicity modes and to focus on 
the features crucial for the Vainshtein mechanism. To introduce with some details this 
mechanism, and illustrate it in action, we start then to discuss it within the framework 
of this decoupling limit, which is much simpler. This framework can also be obtained 
without referring to massive gravity from k-moufiage models [60], generic scalar-tensor 
theories with kinetic self-interactions in the scalar sector, which share many features 
with massive gravity theories, as we will see below. 

5.1.1. Vainshtein mechanism and k-mouflage. Our starting point will be here the 
action 

S = J d^x{ - \h^''Sf,Kp + h^^'X^ij + Mlm^K^L + T^M"]- (67) 

which is a quadratic action for the (normalized) graviton h^y, and describes in addition 
a scalar degree of freedom (f) (which has here the dimension of a mass) with a non linear 
action density -R'atl which is a dimensionless nonlinear function of and its derivatives 
(not necessary only first derivatives) . Note that has a kinetic mixing with h^y through 
the second term on the right hand side above. 

This action can be obtained from the following covariant action 



'>k—mouflage — -''^p I d Xy 



ZIVlp 



Sm[gl (68) 



which defines the k-moufiage family [60] . In the regime of weak gravity we can expand 
(|68ll in perturbations (in terms of the unnormalized graviton h^^). If we do so and 
keep the Einstein-Hilbert term up to h"^, keep only the first order in h in the mixing 
term ~ (pR (this because (f) will be of order of h or of higher order as we will verify 
later), and do not expand Knl (because the Vainshtein mechanism relies precisely on 
the non-linearity of this term), we obtain, after normalizing the spin-2 perturbations (as 
h^u — ^ hfj^i,/Mp), precisely the action ( 1671) . up to a total derivative. 

Coming back to action (|67|1 and redefining the spin-2 mode as h^i, = h^^, — rj^^cj), 
we get 

S = J d^x[-h'^''£^fh^^+^^D4>+M],m.''KNL+^ [t^M" - T$) }.(69) 

^ Note that exact (and proven to be everywhere non singular) solutions appropriate to describe a static 
spherically symmetric space-time on the brane in the DGP model, which would be analogous there to 
the solutions found in [371 liSl IMj , are not known, despite a lot of effort to find them [SSI ISZl [SHI HH] ■ 
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Note that similar to the decouphng hmit of massive gravity, the spin-0 and spin-2 modes 
are decoupled, and a non-minimal scalar-matter coupling appears. The equations of 
motion following from ( 169|) read, 

^,f ^a/5 = ^, (70) 

3n4> + s^ = ^, (71) 

where Sfp = Mpni?5KNL/5(j). The essence of the Vainshtein mechanism can be easily 
seen from the two equations above. First, noticing that when the linear term is dominant 
in the l.h.s. of ( ITTl) . the solution for is of order of /i^y, obtained from (1701) . Hence the 
(normalized) physical metric /i ~ /i + (see Eq. (136|) ) therefore receives 0{1) corrections 
as in a free scalar tensor theory. On the other hand, there is usually a regime where 
the nonlinear term £^ is dominant in fl7T|) . then (f) is subdominant in comparison to h, 
therefore h c:^ h and GR is restored. 

The detailed study of the k-mouflage model for various nonlinear terms K^l in the 
static and spherically symmetric case confirms the arguments above [60]. In this case, 
as well as for what follows, it is convenient to choose the Newtonian gauge for the metric 
written in the spherical coordinates, 

rfs^ = -(1 + ^/Mp)dt^ + (1 - ^/Mp){dr^ + r^dn^). (72) 

The independent components of the linearized Einstein tensor in the gauge (1721) read, 
£oohai3 = ^ (r^$') , £Q{haf3 = 7 (^ — $)'• The linearized GR solution is then 
"^GR = ^GR = —MpTs/r. Typically, the scale where the nonlinear regime switches on 
depends on the form oiEif,. We can parametrize the nonlinear term as E^f, ~ d^~^^^ (p^ / K^, 
with AJ^ = Mpm"'~^. The Vainshtein radius — where the regime changes from linear to 
nonlinear — can be found by comparing the two terms in the l.h.s. of (I7ip and using 
(9 — )■ 1/r, (r being the distance to the source). One finds [60] 



ry = ^{Mprsf-'^/\ (73) 

Note that the galileon model [TJl [IS] is precisely of the k-mouflage type, therefore 
one naturally expects the Vainshtein mechanism to be operating in this model as well. 
Indeed, the study of static spherically symmetric configurations confirms this [151 El] • 

It should be stressed, however, that it is not enough to establish the existence of 
the two regimes: the Vainshtein one, where GR is restored, and the linear one, far 
from the source, where the solution for the linearized theory (e.g. the Fierz-Pauli linear 
massive gravity) is recovered. It should be also checked that there exists an everywhere 
non singular solution that matches the two regimes. This crucial step is lacking in 
the original paper of Vainshtein as was stressed in particular in Ref. [13] and, it is 
only recently that the existence of such a solution was shown [371 ESI El] after some 
previous contradictory claims [621 130] . We will come back below on these results. Let 
us before introduce the Vainshtein mechanism successively in the decoupling limits of 
generic NLFP and dRGT theories. 
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5.1.2. Massive gravity with general potential. In this case, looking only at the 
decoupling limit, one notices that the equations of motion for spin-2 and spin-0 have the 
form identical to those of the decoupling limit of k-mouflage ( 167|) . Indeed, equation ( HB 



coincides with ((TID with £^ ~ 9^0VA^. 
that GR is restored inside ry = {rs/m'^] 
nonlinear differential equation. 



where we have defined /i as 

/i= (p, 



Mprs 



One then expects, from the arguments above, 
^/^. Integrating (HTi) once yields a second order 



(74) 



(75) 



and Q{jx) is a non-linear second-order function of /2, containing terms jljx' , fi'"^, jljl" and 
fi'jji". The function Q{fJi) is in one to one correspondence with the quadratic terms in 
the right hand side of Eq. fBT]) (for details see [SZ])- An integration constant has been 
chosen so as to match the source, resulting in the r.h.s. of (17^ . The quantity /i has here 
the meaning of a gauge function that goes from a unitary gauge where the flat (fiducial) 
metric f^^ has just the Minkowski form in spherical coordinates to a non unitary gauge 
where the only non trivial components of the physical metric are gu and g^- (like in the 
Schwarzschild form) while the fiducial flat metric contains fi [37]. Depending on the 
form of Q{jl) (that is, on the form of the potential term), a global solution matching the 
flat asymptotic behavior and the Vainshtein regime may or may not exist. In particular, 
for the potential f lTSj) there is no such solution. On the other hand, for the potential 
( 1T9|) . such a solution exists. Let us study this case in more detail as an illustrative 
example. For the potential (fT9|) we have 

(76) 
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(77) 



As it was expected, deviations from GR are of order of one in this regime. In particular 
ones notes that, in this regime, \l/ = 2$ which is another way to express the vDVZ 
discontinuity. Note also that we recognize above the flrst terms of an expansion in 
the Newton constant (or equivalently in the Schwarzschild radius of the source). The 
expansion breaks down at the Vainshtein radius where the non linear regimes switches 
on and O (^) ~ 1. In this second regime, we obtain from ( 171|) . 



$ = $ 



Mprs , Mprs ,^ I r 



5/2 



2^/2M, 



prs 



- -7^ 



r r \ry / 9 r 

Therefore GR is restored for r <^ ry as it was anticipated, since we now have \E' = $ 
at dominant order. Moreover, as announced above, the corrections to GR (as well as 
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0) are seen to be oc r^ which is found to be proportional to the square of the mass 
of the graviton m'^. These corrections represent the first terms of an expansion around 
GR solutions in terms of m"^. 

A rigorous study of the Vainshtein mechanism for nonlinear massive gravity in 
the decoupling limit for general potentials can be found in [37]. For some choice of 
the potential in NLFP theory, there are solutions inside ry, which do not feature the 
Vainshtein scaling, however, GR is also restored. This type of non- Vainshtein scaling 
is obtained by expanding (17^ around the zero- mode of the operator Q (176|) . These 
solutions turns out to be unphysical, and it is interesting to compare this to the case of 
dRGT theory, where solutions other than Vainshtein-like also exist. 

5.1.3. dRGT massive gravity As in the case of NLFP with a general potential term, 
the decoupling limit of the dRGT model with the choice /3 = also falls into the k- 
mouflage category. Indeed, the equations of motion for k-mouflage f l70|) - fl7T|l reproduce 
those for the decoupling limit of dRGT model (with /3 = 0) (I65l)-( l66|) with £^ of the 
Galileon type. For spherically symmetric static ansatz, ( 166|1 reads, in terms of fx defined 
as in (1751) 

The last equation, describing the scalar mode of the dRGT model is quite similar to that 
for massive gravity with a general potential ( FMll . with, however, an important difference: 
f l79|) is algebraic, while f l7i|) is a differential equation of the second order. This can be 
traced back to the fact that massive gravity with a general potential possesses an extra 
(BD-ghost) scalar degree of freedom compared to the dRGT model, and that, as we said 
before, the presence of this ghost manifests itself in the higher derivative nature of the 
field equations in the scalar sector, here given by those (as seen in the definition (I75|) ) 
for jl. 

The Vainshtein radius for the dRGT model, indicating where the linear regime for 
this model breaks down, can be read off from ( !73|) and reads 

ry={rslm'y'\ (80) 

The Vainshtein regime is found from ( !79|) assuming the last term in the r.h.s. is dominant 
[63], 

r r \rvj r r \rv J , . 

~ MpTs f r ^ 



^ (2a2)i/3r Vry, 

Note that as in the case of NLFP with general potential, not all a lead to a global 
solution matching flat asymptotic aX r > ry and the Vainshtein regime at r < ry. For 
a > such a solution exist, and for a < it does not llSl 



+ Different conditions for existence and stability however apply when there is cosmological evolution 

mm. 
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When /3 7^ 0, the decouphng hmit of dRGT theory (lFr|) does not exactly fall into 
the k-mouflage category. In this case, indeed, it seems that no local redefinition of 
fields can decouple spin-0 and spin-2 modes [20]. However, it is not difficult to refine 
the argument we used for the k-moufiage model. Indeed, anticipating the Vainshtein 
regime, we neglect the second term on the l.h.s. of ( l65ll and all the terms containing 
only scalar in f l66|) . since we expect h ^ (p. To cancel the r.h.s. of fl66|) . one should 

assume (f) ~ r'^Mpm?/ y jS, so that ^ $, hence GR is restored. More precisely, taking 
\l/ = $ = ^GR and substituting this into f l66|) and neglecting all the terms on r.h.s. but 
the mixing term, we find for the the spherically symmetric ansatz, 

Mprs MpTs ,^ / r \ ^ ~ Mpm"^ , , , 

r r \rvj n 

and ry is defined as in Eq. ( IHOl) . 

For the spherically symmetric static ansatz it is in fact possible to combine (16T!) 
and (162|) to obtain a single quintic algebraic equation, say, on function jl [511 166] . 

2^+2Ai^ +l^Y+3J2Al-96Al = ^^l^'-4Afj' ^^'^ 

and the potentials of the metric are expressed in terms of /i, as 

Note again that Eq. ( 183|) is similar to the equation found for general NLFP theory in 
decoupling limit, ( 171|) . ( 17^ . with the same difference as above (and the same reason 
behind): it is algebraic and not differential. It is worthwhile to mention that since the 
equation for jl is in general of order 5, there exist several branches of solutions, and 
depending on the parameters of the quintic equation, some of them yield real solutions. 
Moreover, as it was found in [51], there is a branch for which inside the Vainshtein radius 
the solution does not approach GR; in fact gravity becomes weaker when approaching 
the source. It turns out, however, that these solutions do not match a fiat asymptotic 
behaviour [67] . 

5.2. The Vainshtein mechanism in complete Non Linear Fierz-Pauli theories: 
recovering General Relativity away from the decoupling limit 

The discussion of the previous section suggest that in untruncated massive gravity 
theories (and a class of scalar tensor theories with non-linear kinetic self- interaction), the 
Vainshtein mechanism works. The results exposed so far for massive gravity, however, 
have been obtained in an approximation scheme: the decoupling limit, in which most 
of the interactions were neglected and it is possible to separate different helicities. One 
could hence ask how these results extend to the complete theory, i.e. when one does not 
truncate to the leading non linear interactions among these helicities. Results, mainly 
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based on numerical works, have been obtained addressing this question as we are going 
to explain here. 

Static, spherically symmetric, and asymptotically flat numerical solutions of NLFP 
with a source were first studied in [30] which found only singular solutions and hence 
concluded that the Vainshtein mechanism was not correct. This was later re-examined 
for the same theories, using more sophisticated numerical methods, in [53l[5l], where, 
for the first time, solutions featuring the Vainshtein recovery were found. The solutions 
show a recovery of the Schwarzschild solution of GR via the Vainshtein mechanism for 
the potential of the form (fTOjl . The ansatz for the physical and the fiducial metrics reads 
respectively, 

f,,dx^dx' = -dt" +(l- U^\ ' e-'^(^)rfr2 + e-^^'\^dn^ , ^^^^ 

so that the physical metric has a form which is convenient for comparison with the GR 
Schwarzschild solution, while the metric f^y has a non-canonical form, even though it 
parametrizes a Minkowski space-time. The function /x appearing above is related to the 
previously introduced fi of the DL as jl = m?Mpix. This ansatz, where both metrics are 
diagonal, is not the most general one for a static and spherically symmetric case (see 
e.g. [68]), however it is the one used to describe the Vainshtein mechanism. For this 
ansatz, the modified Einstein equations and the Bianchi identity (125!) form a set of three 
independent equations. Supplied with the conservation equation for matter source (for 
simplicity, a smoothly distributed source is taken to be described by a perfect fiuid with 
energy density p and pressure P), a set of four quasilinear ODEs for the four independent 
functions A(r), u{r), fi{r), P{r), was obtained, subject to boundary conditions at infinity 
and at the origin. By solving numerically the field equations, using both shooting and 
relaxation methods as well as analytic insights, everywhere non singular asymptotically 
fiat solutions {z/, A,/i, P} were found in Refs. [53l[5l|- These solutions were found there 
to exist for a large range of parameters of the theory and of the source, but such that the 
source is not too compact (i.e. r© > Sr^). They are non singular and asymptotically fiat 
and have the right boundary conditions at the origin r = 0. They feature a recovery of 
GR at distances r -^ ry, where is was also checked numerically that the first correction 
to the GR behavior agree with the analytic expressions obtained by expanding around 
the Schwarzschild solution. They are also well approximated by the solutions obtained 
in the decoupling limit in the expected range of distances rs "^ r '^ m~^, including in 
the linear regime at ry <^ r ^ m~^ where deviations from GR (due to the presence 
of the scalar polarization also responsible for the vDVZ discontinuity) are found. At 
distances above m~^ the solutions feature the exponential Yukawa falloff, which again 
agrees with the the analytic results found by linearization. A crucial feature (which was 
missed in Ref. [30]) of the asymptotic solutions at infinity is that they are not uniquely 
defined by standard perturbation theory (i.e. expanding the solution in the Newton 
constant, as started in fl77j) ). but have non-perturbative hairs allowing to match sources 
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at smaller distances 

Besides the potential f lT9|) other potentials were also studied [M], in particular 
those which were shown in pTj to possess only a non- Vainshtein like scaling at r < ry 
in the decoupling limit. According to the numerical investigations of [51], these solutions 
found in the decoupling limit do not seem to continue into non singular solutions of the 
full field equations. A new limit was also introduced, the weak-field approximation, 
which captures all the salient feature of the solution, including the Yukawa decay and 
the Vainshtein crossover. Although the theories that has been studied in Refs. [531 El] 
suffer from pathologies, in particular the Boulware-Deser ghost problem, these references 
were the first to exhibit an actual and complete proof that the Vainshtein mechanism 
can actually work in a massive gravity theory. These results were later confirmed by an 
independent numerical integration by Volkov p9] and analogous result were obtained in 
the later discovered dRGT model by this author [70] as well as by another group [71] 
which used the shooting method to exhibit solutions featuring the Vainshtein behavior 
at small radii and the linearized Fierz-Pauli behavior at large radii. Another numerical 
investigation of spherically symmetric solutions in dRGT massive gravity was performed 
in [72], where, however, only the case with a very large graviton mass, rsm ~ 1, has 
been studied. A numerical study of the cubic galileon model coupled to gravity [73] 
has shown that the Vainshtein mechanism works in this model, also when taking into 
account full non-linear gravity. 

In the bi-gravity extension of dRGT model |74], numerical investigation of static 
spherically symmetric solutions has been carried out using in particular a multiple 
shooting method [70] (see also Volkov's contribution in this volume). Asymptotically 
flat solutions which exhibit the Vainshtein mechanism of recovery of General Relativity 
at finite distances, in the presence of a matter source have been found. 

Note however that none of the found numerical solutions, either in a generic NLFP, 
or in dRGT model or its bimetric extension, have sources with high compactness. 
Indeed, when one increases the density of the object, the numerics becomes unstable 
and singularities are found to appear [HU [69]. It is still not clear if those singularities 
are physical or if they can be attributed to numerical artifacts. In the first case, it would 
indicate that the solutions cease to exist even before reaching a black hole size. On the 
other hand, one can show that the standard Vainshtein mechanism does not work for 
black holes. Indeed, Ref. [68] shows using geometric arguments, that a geometry such 
as the one given by (1851) where the physical metric g would contain a black hole horizon, 
must necessarily stop at or before the horizon. This results does not rely on the field 
equations and hence applies to a variety of cases including generic NLFP theory, dRGT 
theory and its bimetric extensions. This raises the interesting issue of the end point of 
gravitational collapse in these theories, given in particular that black holes solutions with 
metrics that are not both diagonal (and hence with no standard Vainshtein mechanism) 
are known to exist there. 
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5.3. Time dependent and other situations with Vainshtein screening 

First, following our discussion about static cases, it is worth stressing that although the 
Vainshtein mechanism successfully operates in these cases, it might not be enough to 
pass standard gravity tests in the solar system. Indeed, in a generic shift-symmetric 
k-mouflage model, the Vainshtein mechanism does not appear to screen sufficiently the 
induced time variation of the Newton constant, whenever the scalar field has some 
cosmological dynamics setting the boundary conditions at infinity (this result, first 
obtained in [75], has been confirmed by explicit study of the Horndeski theory [76]). 
Since the variation of the Newton constant is stirringly constrained by the Solar system 
experiments, this puts severe constraints on the direct matter-scalar coupling. Another 
related aspect is that non trivial Galileon profiles have been found to appear even around 
black holes due to a cosmological time evolution of the scalar [77] . 

Another important time-dependent situation is that relevant for gravitational 
waves. Because massive gravity has at least five propagating degrees of freedom, the 
gravitational radiation from binary pulsars should occur at higher rate. One, however, 
could also expect that the Vainshtein screening dump at least the radiation from 
the scalar degree of freedom. As an example of a system possessing the Vainshtein 
mechanism, gravitational radiation from binary pulsars has been studied in a Galileon 
model in Refs. [7H1 [79]. In case of the cubic galileon it was found that the radiation 
is indeed suppressed, although not as much as the static fifth force. This is because, 
firstly, monopole and dipole radiation were found to exist, and secondly, because the 
suppression factor is larger than one could naively guess. Indeed, one could have 
expected that because the static scalar field inside the Vainshtein radius is suppressed 
by a factor (r/ry)^/!^], the radiation is suppressed by the same factor, in comparison 
to the Brans-Dicke theory. However, due to the presence of another scale in the 
problem, the orbital period flp, the suppression turns out to be milder, namely by factor 
{Q,prv)~ [78]. In case of general galileon the study of radiation becomes more subtle, 
since the perturbation scheme breaks down for physically interesting situations [79] . For 
the bi-gravity extension of dRGT model, a recent discussion points out the possibility 
to detect a graviton mass via observations by gravitational wave detectors [80J. A 
possibility to detect gravitational waves from massive gravity by means of the stochastic 
gravitational wave observations has also been studied in [8Tj . 

Other situations of great interest, where the Vainshtein mechanism has to be taken 
into account, are cosmology and many body problems (see in particular Refs [82]). 
However, we will not discuss this here and refer the reader to the contributions of 
Volkov and Koyama in this volume. 

* To get the static spherically symmetric solution in the case of the cubic galileon, one should take 
Eq. ([79]) without the cubic term. Then instead of {r/ry) and (r/ry)'^ factors, in Eq. I|8ip. one obtains 
a (r/ry)'^/^ suppression of the scalar field. 
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6. Conclusions 

The Vainshtein mechanism has been the subject of many recent investigations and 
became a central piece in many theories of modified gravity. Although proposed more 
than 40 years ago, it is only very recently that the original idea of Vainshtein was 
shown explicitly to be correct. Several questions remain however open. Among the 
most pressing ones evocated above, we could mention the nature of the end point of 
gravitational collapse of a star or the way the Vainshtein mechanism can hide sufficiently 
time variations of the Newton's constant in the solar system. More generally, very little 
is known about time dependent cases. To conclude, we would like also to stress that 
the Vainshtein mechanism relies crucially, as it should appear clear from this paper, 
on being able to take into account non linear strong self-interactions, which are in 
general non renormalizable, in some sector of the considered theory. This calls for 
an understanding of UV completion of these theory that is currently cruelly lacking. 
Various possibilities have been put forward concerning this questions in the various 
theories concerned including the DGP model [831 EH ES], but this definitely deserves 
more investigations if one wants to put these theories on a firm footing. 
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